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>< : ABSTRACT 

A set of different conformal solutions corresponding to a constant flux of 
squared vorticity is considered. Requiring constant fluxes of all inviscid vorticity 
invariants (higher powers of the vorticity), we come to the conclusion that the gen- 
eral turbulence spectrum should be given by Kraichnan's expression E(k) oc A; -3 . 
This spectrum, in particular, can be obtained as a limit of some subsequences of 
the conformal solutions. 
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The problem of small-scale spectrum of two-dimensional turbulence is a pecu- 
liar problem among the variety of turbulent systems. The point is that dimensional 
considerations do not give a steady spectrum that corresponds to the enstrophy 
(squared vorticity) cascade. The spectrum obtained from dimensional analysis is 
E{k) oc A;~ 3 [1] which yields a logarithmic infrared divergence after substitution 
into the equations for the correlation functions. Since any cascade picture assumes 
turbulence locality the divergence makes this spectrum rather suspicious. Kraich- 
nan's attempt to save the spectrum from nonlocality by introducing the slow factor 
ln _1//3 k attains convergence only in the first order of perturbation theory [2] while 
the next orders reveal divergencies with higher powers of the logarithm: In etc. 
The fact that the powers of the logarithm increase with the order of perturbation 
theory suggests that a substantial renormalization of the index occurs. No suc- 
cessful attempts to work out the divergencies or to show that they are cancelled 
are known to us. The existence of alternative predictions for the steady spectrum 
E{k) oc k~ 4 by Saffman [3] and E(k) oc Ar 11 / 3 by Moffatt [4], show that this is 
still an open problem. 

A fairly new approach to the problem has recently been introduced by Polyakov 
[5]. He suggested to borrow a set of correlation functions from conformal field 
theory to satisfy a chain of equations following from Euler's equation. Note that 
the conformal invariance of turbulence should be considered as a pure conjecture. 
It is convenient to start from the Navier-Stokes equation written for the vorticity 
field u;(x, t) 

duo dtp dAvp ^ 

dt dx a dx/3 

Here if) is a stream function giving the velocity field as follows: v a = e a pdpip. Here 
and below we use a shorthand notation d/dx a = d a . 

Our aim is to find a stationary set of equal time correlation functions 

J n (xi, . . . , x n ) = (w(xi, t) . . . u(Xn, *)) • 

The brackets denote an average with some time independent probability distribu- 
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tion: 

i:(-(x l! t)...^^...c,(x, l) t)\=0. (2) 
p=l ' ' 

Such a stationary set is expected to exist in the inertial interval of scales, i.e., for 
distances that are much less than the scale of an external pump (or the scale of 
an initial distribution for free decay) and much larger than the viscous scale a. It 
is possible, then, to neglect viscosity in (1) using instead a careful point splitting 
procedure [5] for the nonlinear term 

g djjjx) dAiP(x) _ c <9?/;(x + a/2) <9A^(x - a/2) 

Q ^ dx a Oxp Q/3 dx a dxp 

where "lim" implies angle averaging. To calculate different-point pair correlators 
like ip(x. + a)^(x — a), the fusion rule of the type 

m m = M + • • • (4) 

should be used. Here we follow the notations of Ref.5 so that [ip] means the 
conformal class of if), i.e. itself together with the operators L_ m . . . L^ rik ip, L_ n 
being Virasoro generators [6]. Both ip and are presumed to be taken from a set 
(primary fields) of some conformal field theory (in this paper we will consider the 
so called minimal models [6]). The primary field (f> provides the main contribution 
in the operator product expansion (OPE) (4) in the small-scale region, i.e. it has 
the smallest conformal dimension. The important thing is that the scaling indices 
(dimensions) of the fields are not additive so generally ^ 2A^. The energy 
density in the wave number space is expressed via |^| 2 and is 

E(k) = 2ixke(k) oc k 4A * +1 , (5) 

while the enstrophy density is H(k) = k 2 E(k). Here e(k) is the energy density in 
k-space. 
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To choose an appropriate solution from the wealth of conformal solutions, 
one should impose some additional conditions that follow from the symmetries 
or conservation laws specific for the problem in question. According to Fjortoft's 
theorem (see e.g. [7]), the vorticity is the relevant quantity in the problem of 
small-scale turbulence (while the energy flux determines large-scale turbulence). 
Following Kraichnan [1] who developed a simple and efficient (though uncontrol- 
lable, of course) closure in terms of double correlation functions, the enstrophy 



H 2 



= J u 2 (x)d 2 x , (6) 



which is a motion integral of Euler's equation is usually taken into account. A 
steady turbulence spectrum in the small-scale region should provide for a constant 
enstrophy flux over the scales which yields [5] 

^^±^^))ocr° = const. (7) 

Puting uj = Aip and 

for the time derivative, Polyakov obtained [5] 

(A^ + 2) + (A^ + l) = 0. (8) 

As one can see, the enstrophy flux is expressed through the triple correlation func- 
tion which can be expressed by the fusion rule (4) in terms of the double correlation 
function. 

Equation (8) can be obtained also by requiring the rate of the enstrophy dis- 
sipation to remain constant while the viscosity v goes to zero: 

l/a 

^ = v ] k 2 H k dk oc ua- 6 ^ oc z/ (3+A,+A,)/(A,-A,) (Q) 
The last estimate was given by using the expression for the viscous scale v oc 
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a * * that follows from the comparison of the nonlinear and the viscous terms 
in the Navier-Stokes equation. 

Kraichnan's dimensional approach would correspond to additive dimensions 
A^ = 2A^, giving thus A^ = — 1 and E{k) oc k~ s . However, we have arguments 
that suggest that this in not a conformal solution from the set of minimal models 
(see Appendix). 

To ensure that the conformal set of correlators is a steady solution, Polyakov 
imposed an extra condition requiring that rhs of (4) (which determines the time 
derivatives of the correlators) vanishes in the ultraviolet limit. Since 

*P(ziMz 2 ) = (*i - z 2 ) A ^ 2A ^(z 2 ) + ... 

as z\ — ■> z% the following inequality arises: 

A > 2A^ , (10) 

thus, using (8), < — 1. 

Conditions (8) and (10) by no means determine a single solution. The minimal 
model (2,21) presented by Polyakov [5] is nothing but the first example from an 
infinite family. The curious reader can find the first few hundred solutions in the 
Appendix. Polyakov's solution corresponds to the minimal number of primary 
fields (in this case 10). One could, in particular, find the minimal model (5,72) 
that gives A^ = An 2 5) = —7/6 and E(k) oc A; -11 / 3 as in Moffatt's spectrum. 

Usually, when speaking about a turbulent solution carrying a constant flux, 
one should check that two conditions are satisfied: i) The solution should be lo- 
cal in A;-space which means that distant scales do not interact substantially; this 
should be provided by the convergence of the integral determining the flux in k- 
space; ii) The constant that arise in this (converging) integral should be nonzero 
and have the correct sign to satisfy the boundary conditions in /c-space, i.e., the 
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pumping and damping. Any solution of (8) and (10) violates both of these condi- 
tions. Inequality (10) means that Ay, < —1 i.e. the spectrum (5) is steeper than 
Kraichnan's one (which yields an infrared logarithmic divergence) so that a power 
infrared divergence arises for any solution. (The correct power of the divergence 
can be obtained by using Lagrangian or quasi-Lagrangian variables which eliminate 
sweeping of small scales by larger ones [8,9]). For minimal models the three-point 
function {ipipip) ~ {(pip) is equal to zero since primary fields with different dimen- 
sions (Ay, 7^ A^) are orthogonal. The cases where ip appears in the operator 
product expansion of [ip}[ip] (like those with Ay, = A^ = — 3/2 which one can find 
in the Appendix) break parity and therefore should be excluded. The conformal 
solutions in question are thus fluxless. The second difficulty (zero flux) remedies 
to some extent the first difficulty (nonlocality) since one should not require the 
convergence of the integrals that are identically zero. Physically this corresponds 
to the fact that the spectrum of a system in thermodynamic equilibrium should 
not be local, unlike the cascade spectrum one. A fluxless spectrum corresponds to 
an equilibrium case. 

But, since we are discussing a nonequilibrium situation, there still remains the 
question, how the spectrum carries nonzero flux from the pumping region to the 
viscous region of scales. Polyakov's suggestion is that small deviations from the 
power law due to infrared cut-off (pumping scale) could provide nonzero values for 
the enstrophy flux. (It is unclear whether such a flux is local in /c-space or it is 
a noncascade nonlocal solution.) The physical correlators are thus assumed to be 
close to the equilibrium ones in the inertial interval of scales. This is similar to 
what happens in two-dimensional optical turbulence described by the Nonlinear 
Schrodinger Equation 

im t _|_ + T|^| 2 = 0. 

For wave turbulence, the small-scale turbulent spectrum carrying constant energy 
flux is as follows: 

e(Jfe) oc k a - m - d , 
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where d is the space dimension and a and m are the scaling indices of the Hamil- 
tonian coefficients (i.e. the frequency and the four-wave interaction coefficient 
respectively) [10]. For the NSE, a — d — 2 and m = so that the turbulent 
spectrum is e(k) = const which coincides with the equilibrium equipartition. This 
spectrum (which is an exact steady solution) is fluxless too. Numerical simulations 
of the NSE show the nonequilibrium spectrum to be close to e ^ const, an analyt- 
ical attempt to introduce a slow logarithmic factor causes some doubts [11]. The 
same coincidence of the equilibrium and turbulent spectra takes place for common 
turbulence of Langmuir and ion sound waves in plasma, the spectra carrying fluxes 
acquire logarithmic factors in this case [10,12]. Polyakov suggested to distort the 
spectra by analytical (in x-space) contributions. In our opinion, the main differ- 
ence of his approach to hydrodynamic turbulence from the above picture of wave 
turbulence, namely, the assumption that the turbulent spectrum of 2d hydrody- 
namics should be close to an equilibrium one, is not based on solid ground. One 
would like to see the degeneracy that prescribes the coincidence of the turbulence 
spectrum with an equilibrium one. 

And what may be more important, Polyakov's considerations do not take into 
account the presence of an infinite set of motion integrals 

H n = J U n (x)d 2 X . (11) 

The conservation of H n follows directly from the fact that the equation 

— -e — — — = (12) 

dt a ^ dx a dxp 5uj 

conserves the integral J F(ui) dxdy, where F is an arbitrary function and the Hamil- 
tonian 7i is an arbitrary functional of oj [not only 7i — J ipcu dxdy, which gives the 
lhs of (1)]. As one can see, even an infinite number of motion integrals does not 
fix the system but only its class. 
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Most authors feign indifference to the existence of the infinite number of motion 
integrals in 2d turbulence. Some (weak) arguments that only quadratic integrals 
(i.e. energy and squared vorticity) should be taken into account while considering 
thermodynamic equilibrium were given by Kraichnan [1,12]. However, an arbitrary 
turbulent pump generally produces a nonzero input of all integrals H n . The theory 
should describe the fate of these integrals. Estimating the rate of the viscous 
dissipation of H n (n > 2) similarly to (9), one gets 



Formula (14) again demonstrates that < —1 while — A^ < 0, otherwise 
the dissipation rate goes to infinity while viscosity goes to zero. For Polyakov's 
solution with A^ < — 1, the integrals H n are not dissipated in the inviscid limit 
when n > 2. What is the fate of these integrals if they are injected? Note that it is 
impossible to have a local cascade of H2 and a nonlocal transfer of other vorticity 
integrals. 

One can require a constant dissipation rate of any vorticity integral. According 
to (13) and (14) it means that A^ = —1 and we are coming back to Kraichnan's 
spectrum. Indeed, this corresponds to ip(r) oc |r| 2 so that the vorticity u has 
zero scaling index (maybe logarithmic). All powers of the vorticity can thus have 
constant fluxes in /c-space simultaneously. Actually, Kraichnan's spectrum equally 
satisfies all conservation laws. 

Returning to the above set of conformal solutions, one can find subsequences 
(see Appendix) that give A^ which approaches the value —1, while the number 
of primary fields in the model increases. The closer A^, is to —1, the better the 



1/0 




(13) 



We rewrite it by the help of (8) as follows: 



l/a 




(14) 
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given conformal solution satisfies the vorticity conservation laws. One can suggest 
to get Kraichnan's spectrum as a limit of such subsequences, thus considering it as 
some weak solution. In hydrodynamics, it is physically intuitive, to account for an 
infinite rather than a finite number of primary fields (of topologically different field 
configurations). Still, some other spectra (including those of Saffman and Moffatt) 
could be realized for special initial conditions or as intermediate time asymptotics. 
Kraichnan's spectrum seems to correspond to the most general conditions (both 
initial conditions in time and boundary conditions in /c-space). This spectrum 
seems to be the most ergodic (e.g. isovorticity lines have the fractal dimension 2 
so that they can fill the entire space [13]). It is quite natural that Kraichnan's 
spectrum corresponds to an infinite number of primary fields. How this spectrum 
can be distorted by a slow factor to avoid the logarithmic divergence is, in our 
opinion, still an open problem. 

Acknowledgements: Discussions with P. Wiegmann, A.Finkelstein, G.Schutz, R.Plesser, 
A.Schwimmer and Y.Levinson are gratefully acknowledged. 

APPENDIX 

In this appendix we will summarize all the conditions set by Polyakov on the 
solutions of a turbulence problem. Subsequently we apply these conditions to the 
set of minimal models and present a list of a few hundred solutions. We describe 
the algorithm we used to compute these solutions and give arguments why we think 
A^, = — 1 is not a minimal model solution. Finally we present a few solutions in a 
sequence that tends to A^, = — 1. 

The conditions are: 

1. A^ + A = -3 

2. A^ < -1 

3. A^ must be the operator with the smallest dimension in the OPE of [tp] and 
[# 
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The minimal models are characterized by two positive co-prime integers (p, q) . 
For each minimal model (p, q) there is a set of primary fields parame- 

terized by two integers (n, m), 1 < n < p — 1, 1 < m < g — 1. The spectrum of 
conformal dimensions is given by 

{nq-mp) 2 - (g-p) 2 
Vm) = 4^ • (A- 1 ) 

From this formula we see that A( njm ) = A( p _ ng _ m ) and that they both correspond 
to the same primary field. The scaling dimension of the primary field (n, m) is 
2A( n m ). From the first condition, using (A.l), we get | in terms of n^, m^, n^, 



p _ n^m^p + n^m^ - 8 + k 
q m 2 + m| - 2 



(A.2) 



where k is an integer defined by 



k 2 = {n^m^ + n^m^ - 8) 2 - (n 2 + - 2) (raj + raj - 2). (A.3) 

This defines both p and q since they are co-prime. From the third condition and 
the selection rules for the OPE we get that ra^ra^ are odd numbers and satisfy 
1 < n<p < 2n^ — 1, 1 < < 2m^ — 1. From the demand that be the lowest 
possible we get {n^q — m^p) m 0, this sets ra^ to be the nearest odd to 

n^ra^ - 8 ± w2m 2 - Wn^m^ + 2raJ + 60 
n V -g— . (A.4) 

'V 

An algorithm to calculate these turbulent solutions can be as follows. Given 
n^, m^, one can calculate ra^ from (A.4) and if k (eq. (A.3) ) is an integer one can 
olso find p and q. Using this algorithm we got a set of solutions which are listed in 
the tables below. (This is a partial list, there is an infinite number of solutions.) 
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We also want to check if there exists a minimal model for which Aw, = — 1 or 
A^ = —2. We set |n^g — vn^p\ to its minimal value and since p, q are co-prime 
there exist n^, such that n^q — m^p = 1 and we get, using (A.l), —2 = — l P q 
or 

p 2 - Ittpq + q 2 - 1 = 0. (A.5) 
for large p, q one can neglect 1 and solve the quadratic equation to get 

Q ± = 5±V24, (A.6) 

the two solutions correspond to the symmetry between p and q. (This solution is 
irrational so in practice one can take p, q co-prime such that | is close to this value, 
one can get any accuracy by taking large p, q). From A^, = — 1 and (A.6) we have 
in terms of as the nearest integer to 

™>ih — n^- + 2, /-. (A.7) 
p y p 

The sequence defined by 

a/c+i = 10a fc - a/c-i a = 0, a\ = 1 (A.8) 

serves as a solution to (A.5) taking p = a^, q = a-k+i for any A: = 2,3,... the fc-th 
term is given by 

a fc = ( 5 + ^-l 5 -^ fc (A.9) 
2V24 V ; 



Demanding that this (p,q) model has Aw, = —1 leads to the condition that 



y/Apq + 1 is an integer or, in terms of this sequence, that 

(5 + y^) 2 ^ 1 + (5 - V2A) 2k - 1 + 14 
24 

is an integer squared, which is unlikely. 



(A.10) 
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In the list below one can find the models with (p, q) close to (A. 6) and {n$, m^) 
satisfying (A.7): (13,129), (39,389), (109,1082), (232,2295) and (69,686) with 

_A , - 561 391 59001 17873 7905 rp( . n ppt; V pl v 
^-V ~~ 559' 389' 58969' 17748' 7889 ^pecwveiy. 
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